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. Reconsidering the variational procedure for uniaxial systems modeled by continuous free energy 

functionals, we derive new general conditions for thermodynamic extrema. The utility of these con- 
ditions is briefly illustrated on the models for the classes I and II of incommensurate-commensurate 
systems. 
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Numerous materials which are under intense investigations in the contemporary condensed matter physics are 
thermodynamically one-dimensional. The well-known examples are various uniaxial materials with incommensurate 
, and commensurate orderings and quasi one-dimensional conductors with charge or spin density wave instabilities 
[Q. Order parameters for such systems are generally multicomponent, u = (ui,U2, ■•■•,MAf), and depend on a single 
spatial variable x. The principal task is then to find thermodynamically stable configurations Uc{x), those which 
minimize the free energy functional T. Since the latter is the one-dimensional integral, it is tempting to treat this 
', variational problem as an equivalent to the standard classical mechanical one with the roles of time variable, 
^ ' vectors in the N-dimensional mechanical configuration space, action functional and Lagrangian attributed to x, u, J- 
O \ and / respectively, the latter being the free energy density. 

In the present Letter we do not follow this widely accepted attitude, but start from two obvious, yet substantial, 
\ differences between these two variational schemes. The first one is present in the very extremalization procedure. 
' In contrast to the classical mechanical trajectories, the realizable solutions of the Euler-Lagrange (EL) equations 
for thermodynamic problems follow after an additional extremalization with respect to the initial (or boundary) 
' conditions. The second difference concerns the content of the free energy densities. In the most interesting models 
' for incommensurate-commensurate (IC) systems, including the basic ones, they contain either terms linear in the first 
derivatives u — (u'j^, U2, u^), or terms with higher derivatives u^^^ = (j > 1) (or both), in contrast to the 

I standard mechanical Lagrangians which do not contain analogues of such terms. 

' Starting from the first observation, we reformulate the procedure of thermodynamic extremalization, and derive, 

under assumptions specified below, the following necessary conditions for any thermodynamic extremum Uc: 
' Condition A; 



c : ^ 
o 

o 



dx = 0, (1) 



where n is the order of highest derivative of u present in the free energy functional, and L is the length of the system 
^ ■ taken in the thermodynamic limit L ^ 00. In particular, for free energy densities which do not depend explicitly on 
\ X the condition (|l|) reduces to the simple equality 

S ■ + H = 0, (2) 

where is the averaged value of free energy and Ji is the integral constant which has the meaning of Hamiltonian 
in the equivalent classical mechanical problem (but does not have a direct physical meaning in the thermodynamic 
counterpart) . 
Conditions B; 

L 



\l y.-'S^dx = , (3) 



c.a 
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where a = 1, N . 

The ensuing discussion will show that in the case of thermodynamic functionals of the standard " mechanical" form 
the conditions (||) and (||) are of almost trivial meaning. They however have far-reaching implications just in IC 
models, for which, as was already pointed out, free energy densities depend in more complex ways on derivatives u*-"^. 
These conditions also appear to be a powerful tool in the numerical determination of phase diagrams, particularly for 
systems with nonintegrable free energy functionals. 

In order to derive the Conditions A and B we start from the general expression for the free energy functional 



u(j;),u'(a;),u"(a;),...,u(")(a;);a;j dx, (4) 

where / is an analytical function of its arguments, bounded from below. Each thermodynamic extremum Uc{x) of 
this functional has to obey the variational condition 6J-'{{uc}) = 0, equivalent to the Hamilton variational principle 
in classical mechanics. This necessary condition leads to the EL equations 



(-1)^^7:^ = (« = i,...,iv), (5) 



^' dx^ BnM) 



equivalent to the Lagrange equations in classical mechanics. The solutions of the EL equations (||) form a set 
{u(x;./l)} which generally depends on 2nN continuous parameters (ai, a2„7v) = -4. There is a freedom in the 
definition of the parameters A, the most usual choices being initial conditions [u(a;o), u'(xo), u'^"^^^ (^^o)] where xq is 
an arbitrary initial spatial position, and boundary conditions [u(a;i), u'(a;i), ...u^""^^ (xi); u(a;2), u'(a;2), ...u^""^^ (X2)] 
where xi and X2 are arbitrary end points. In classical mechanics these two choices correspond to the Newton and 
the Hamilton (variational) axiomatizations, respectively. Thermodynamic extrema, including thermodynamically 
stable configurations for which S'^T > 0, are those members of the set {u(a;;„4)} which extremalize the free energy 
^({u;^)}) as a function of the parameters A. This additional property completes, together with the EL equations 
(^, the sufficient condition for thermodynamic extrema. In particular, a configuration which fulfills the conditions 
6J- = and 6^ J-' > is thermodynamically stable only if it is also a minimum in the set {u(x; A)}. 

The dependence of J-{{u; A)}) on the parameters A is generally intricate. It may be at least partly nonanalytic, 
as is usually the case for the functionals (^) with nonintegrable EL equations Q, and in particular for those with 
free energy densities / which are explicitly x-dependent. Thus, there is no efficient general way to extract local 
extrema of J- from the set {u(a;;„4)}. However, we can now conveniently reformulate the above proposition that the 
thermodynamic extrema follow from the succession of the first order variation and the extremalization with respect 
to the parameters A, into an equivalent, and again sufficient, requirement that the solutions of the EL equations are 
thermodynamic extrema if they are local extrema in the set {u(a;)} of all configurations allowed by the functional 
(H). By this enlargement of the set within which we are looking for the local extrema Uc{x), we get a freedom to 
choose arbitrarily (and suitably) the parameters with respect to which the set {u(x)} is analytic and corresponding 
extremalizations reduce to simple differentiations. This freedom will be here partly exploited, by making two choices 
of continuous parameters which will lead to the Conditions A and B. 

The first continuous parameter is introduced in the following way. Let us take one thermodynamic extremum, 
Uc(x), and define a set of functions {u(x; q)} by 

u{x;q) =Uc{qx). (6) 
The free energy functional (^ for this set becomes a function of q given by 

J^{{u{x;q)})^F{q)^^ \u,{z), qu'^{z), q'u'^{z), q^ui^Hz); q-^ z] dz, (7) 







with z = qx and u'J\z) = dui{z)/dz^. The requirement that Uc{x) is an extremum in the set {u{x;q)} is expressed 

by 

[dFiq)/dq]g=i = 0, (8) 

provided F{q) is a smooth function of q for q ^ 1. Let us also take the thermodynamic limit L 00 and assume 
that F{q) then does not depend on L [up to the corrections of the order 0{1/L)]. Under these assumptions, which 
will be critically examined later on, F(q) may depend on q only through the density / in Eqj^. The latter is an 



2 



analytic function of q since it is analytic with respect to u', ...,u("^ by assumption. The derivative dF{q)/dq is then 
well defined and the requirement (||), applied onto the function gives the Condition A, Eqj^. 

The further simplification takes place for the functionals (^) in which the free energy density does not depend 
explicitly on x. Then, like in classical mechanics, there exists an integral constant (Hamiltonian) , 
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for each solution of the EL equations (||). Using the obvious identity H ~ H dx, and the identity 
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u'^'^+^^g) dx 



(10) 



which follows after I successive partial integrations of the left-hand side, one reduces the expression (Bl) to 
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(11) 



Here g in Eq.|l^ is identified with df /du"a from Eq.||, and Fc = J-{{vic})- The second term on the right hand side 
in Eq.^ is negligible in the limit L oo, provided Uc(x) and its derivatives are finite. All thcrmodynamically 
stable extrema have this property since / is bounded from below. The third term vanishes for each thermodynamic 
extremum due to the condition (0). The expression (11) thus reduces to the Condition j4, EqJ|. 

The equality (|) is the consequence of the invariance of the functional (^) with respect to translations in x, and of 
its non-invariance with respect to the changes of z-scale. Like in classical mechanics, the former invariance ensures 
the existence of the integral constant H and the degeneracy of the solutions of EL equations with respect to the choice 
of "initial position" xq. The number of parameters on which the set {u(a::;yl)} explicitly depends is then 2nN — 1. 
Note that for all nontrivial functionals (^) one has > 1 and n > 1, so that 2nN — 1 > 1. In the simplest nontrivial 
case = n = 1 the set A has one parameter, i. e., just H . 

For functionals (^ with an explicit x-dependence of /, the insertion of the EL equations (|^) into the expression (||) 
leads to the relation = —H{L), where H{L) is given by the, now x-dependent, expression (^ at x = L. Since the 
right-hand side in this relation depends on L, it is inconsistent with at least one of two assumptions on the analyticity 
of F{q) specified below Eq.||. We come to the conclusion that whenever the free energy density depends explicitly on 
X, all thermodynamic extrema are isolated nonanalytical points of the corresponding functional (^) with respect to 
changes of x-scale. This fundamental property is the reason why the Condition A does not hold for such functionals. 

Our second choice of continuous parameters from the set {u(x)} is defined by the scaling Ua —> SaUa for any 
1 < a < A^. The steps equivalent to those specified by Eqs.^||can be repeated now for each a for which F{sa) is a 
smooth function. The corresponding conditions 



[dF{s^)/dso 







then reduce to the Conditions B. Performing partial integrations and inserting EL equations 
gets the conditions 



(12) 

into Eq.0, one finally 
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(13) 



which are constraints on the boundary values of the thermodynamic extrema. Note that the boundary ("surface") 
terms are here the leading ones, in contrast to the Condition A in which the analogous terms are only negligible 
0{\/ L) corrections to the finite volume terms of the order 0{L^). Obviously, any periodic solution of EL equations 
satisfies the conditions (|l^). To this end it suffices to take into account corrections of the order 0{1/L) coming from 
the boundary terms in Conditions A and B, in particular a correction which adjusts the period to be a divisor of 
L with an integer ratio. No analogous adjustment for the quasiperiodic and nonperiodic solutions is apparent. The 
conditions (^) are therefore expected to represent restrictive constraints on these solutions as possible candidates for 
thermodynamic configurations. 
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The extremalization of the thermodynamic functional (Q) with respect to the parameter set A, and its non-invariance 
with respect to the transformations x qx and particular, become short of physical justification when 

transposed to its mechanical counterpart. For free energy densities which have the form of conservative Lagrangians 
one has n = 1, and the derivatives , Zi2, enter only through a positive definite quadratic form ("kinetic 
energy"). The criterion (|^) then singles out only equilibrium points (homogeneous configurations) u = cte as possible 
extrema. For such solutions the condition (||) is trivially fulfilled, since \\.c{qx) and the corresponding free energy F{q) 
does not depend on q. The same is true for the Conditions B which reduce to [waUalo — CI- 
As was announced in the introduction, the utility of the Conditions A and B becomes apparent for the functionals 
(^ which have richer dependences on the derivatives of u and allow for the x-dependent stable configurations. For 
illustrations we take the basic models for the classes I and II of IC systems |l| , defined by Q 

/ = 1(0' -5)2 -1/(0) (14) 

and 1^ 

/ ^ {u" f {u' f + Xu" + ^u" (15) 

respectively. 

The decisive term in the model (^4|) is the Lifshitz invariant 50'. is the phase variable, so that 1^(0 + 27r) — V^(0), 
the simplest choice being the sine-Gordon model with a single Umklapp term, V{(j)) oc cosp0, where p is an integer. 
The problem (|^, |lj) is entirely solvable |||,|^ , since the corresponding EL equation is integrable and the set A has one 
parameter, e. g., H . Here we show how the Condition A enables an elegant derivation and an original interpretation 
of the solution. The Condition A for the functional (^, |l^ reduces to 

2ir n2-K 



2nS = I,= / 0^(0)d0= / ^2[-F, + 52/2 - F(0)]d0. (16) 



The determination of the thermodynamic phase diagram, i. e., of the dependence of Fc on the control parameters 
present in the model ([l^), is thus reduced to the calculation of the integral Ic- The relation (|6|) also states that for 
a thermodynamic extremum the corresponding mechanical action variable is just equal to the Lifshitz parameter 6\ 
The dependence of the period P of the stable configuration on control parameters follows from the known relation 
for mechanical systems with one degree of freedom, P = dIc/dH Finally, the corresponding configuration (f>c{x) 
follows from the quadrature of the EL equation with an already determined value of H . Thus, using the equality (g) 
we avoid a more tedious procedure used in the analyses of the models (^4|) ||J^, namely, the entire integration of the 
EL equation (with free H) followed by the minimization of the free energy F(H) as a function of H. 

Since the transformation ^ already exhausts the freedom in the choice of variational parameters A for the model 
(p^, the Condition B which is now given by [0(0' — S)]q = 0, cannot be an additional constraint, but may only 
reproduce some already derived property of the extremum 0c. This condition states that the configuration 0c has a 
slope 0' = (5 at the points x = 0,P, 2P, NP, where P is a period and is a large (macroscopic) integer. It indeed 
follows independently from the EL equation and the condition (|l6|). 

Various criteria suggest that the model (|l^) is nonintegrable due to the presence of the second derivative of the 
real order parameter u. Very probably H — {u"f — {u' f — 2u'u"' — Xu^ — ^u'^ is the only integral constant among 
three parameters in the set A. The Condition A now reads 

2 



2("c) 



dx = 0. (17) 



Without using this condition, we have minimized numerically the functional ^q ) in the Fourier basis and showed 
that the phase diagram contains an enumerable set of metastable periodic solutions with homogeneous domains 
connected by sinusoidal segments Q . The subsequent check jl^ verifies that all these solutions satisfy the condition 
(p!?!). Furthermore, by using it, one significantly facilitates the numerical calculation of (meta)stable configurations for 
the model (p^. Namely, the search for local minima in the Fourier basis gives, as a rule, continuous families of periodic 
configurations. In order to find the proper thermodynamic configurations within one family it suffices to determine 
zeros of the diagonal quadratic form of Fourier components to which the left-hand side of equation ( pT[ ) reduces. By 
this we directly confirm that the obtained configuration satisfies the EL equation and determine its period. The more 
detailed presentation of this procedure for the model (|l^) and its various extensions is given elsewhere pO| |. 
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Applying the Condition B to the model ( [1^ ) we obtain an additional constraint on the boundary points, 

[u'u" - u{u' + u"')]o = 0, (18) 

which, together with the arguments given below Eq.|l^, reinforces the expectation based on the independent numerical 
analysis ^ t hat all thermodynamic extrema of the problem ( [l5| ) are very probably periodic. Note that by conditions 
( ]l7| ) and (|18|) we have fixed two out of three parameters from the set A for the problem (|l5|). Very probably T is not 
analytic for any choice of the remaining third variational parameter, in close connection with the nonintegrability of 
the EL equation and the corresponding chaotic structure of the portrait in the phase space {u, u' , u" , u'"). 

Having these and other examples in mind, we connect the limitations of the present method with the degree 
of the nonintegrability of a given functional by the following conjecture: larger is the number of missing integral 
constants (in the classical mechanical sense), smaller is the number of analytic conditions for the thermodynamic 
extrema (like those given by Conditions A and B). 

In conclusion, necessary conditions for uniaxial thermodynamic extrema are obtained from the extremalization 
with respect to space and order parameter scales. This procedure proves to be feasible for the free energy densities 
which are not explicitly dependent on the space variable. In particular, we show that in this case the sum of the 
averaged free energy and the integral constant (Hamiltonian) vanishes for each thermodynamic extremum. Besides 
their general significance, the present results will be certainly of practical use in analytical and numerical analyses of 
particular models for uniaxial systems. 

We acknowledge discussions with M. Latkovic. The work is supported by the Ministry of Science and Technology 
of the Republic of Croatia through the project no. 119201. 
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